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Abstract 

In this letter we research on the mean extinction-time which is the expected time when the value of the asset reach zero in the 
classic Black-Sholes model. We resolve the stationary backward Kolmogorov equation and make a direct comparison between 
our predictions and the numerical simulations of stochastic differential equations. 
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Introduction 

The famous Black-Sholes model was the starting point of a new financial industry, and it has been a very important 
pillar of options trading since then. Throughout this letter S (t) denote the value of an asset at time t that may take 
any non-negative value. In the classic model for asset price behavior (see [3], and Chapters 6 and 7 of [7]), S (t) 
follows a geometric Brownian motion satisfying a linear SDE: 

dS (t) = pS (t)dt + oS (t)dW(t), (1) 

where p and a are two non-negative constants, p is called the drift, and a is the volatility that quantifies the risk. 
Typical values for a lie between 0.05 and 0.5, i.e., a volatility between 5% and 50%, and the drift parameter is 
typically between 0.01 and 0.1. 

It is possible to solve exactly (1) by means of Ito's formula with F(S, t) = ln(S ) (see for instance [1, 11, 8, 7]). Its 
solution is 

S (t) = S (0) exp ((3t + aW(t)), 

where [3 = p - a A 2/2, and its expectation is 

IE(S(0) = S(0)<^. 

There exist two different asymptotic behaviors: 

(i) If p < 0, then Prob j lim \S(t)\ = o| = I . 

(ii) If/? > 0, then Prob \ lim |5(f}| = oo j- = L 

Therefore, according to this model, if the volatility is sufficiently large, the asset price will eventually decay to zero 
with probability one, but without specifying how quickly it will do it depending on a. In what follows, we are 
interested in Case (i). 

Let T fg inf{t > 0 : S (t) = 0} be the random variable that indicates the extinction-time, obviously T depends on the 
initial population size and we denote this dependenceby TS(0), whereas the expectation time until extinction, or 
mean extinction- time, is xS (0) = E(TS(0)). This mean extinction- time satisfies the following boundary value 
problem (BVP): 


( 2 ) 

(3) 




cIt 1 

^Ts + 2 

r(0) = 0. 

-f-(M) = 0 , 

ds 


d 2 T 
ds 2 


= - 1 , 


(4) 
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d 2 r dr 


(5) 


t(z) = ri + C 2 e 






supposing that S (t) cannot exceed a value M. Now, we have a linear differential equation of variable coefficients; 
after the change of variable z = In s we obtain 


u T 


dr 




(5) 


with a = o A 2/2 > 0, whose exact solution is 

Tiz) = C] + Cie^ ah - -j-£ 
jS 


or equivalenly 

t(s) = c\ 4- c 2 s-&* — - In s. 

P 

Nonetheless, it is evident that this last function cannot solve the boundary problem (4), because it is always 
unbounded; therefore, (4) has no solution. To avoid this difficulty, we study the time T at which S (T) = > 0, i.e., 

we will consider the mean extinction time for the new stochastic variable S~ (t) = S (t) - . This variable satisfies 

dS(t) = + t-Jrf/ + <r($(t) + 


and its mean extinction- time x( ~s) = E(T) satisfies 


.dr I j ,2 
r(s + € )— + -cr {s + e) 

ds 2 ds- 

t(Q) = 0 , 

= o. 


= - U 


(6) 


supposing that S~ (t) cannot exceed a value M. We apply the change of variable z =ln( S~ + ) to this linear 

differential equation of variable coefficients, in order to obtain again (5). In particular, without lose of generality, 
we take = 1, for which the solution 

m = ^(M + 1 f a [l - a + 1)-*“] - f ln(j + I). (7) 

solution that tends to x( ~s) = ln( ~s + 1)/ 1 [3 1 , as M — » °°, because [3 < 0 and a > 0. 

On the other hand, as we announced in the introduction, we have simulated (1) using the Euler-Maruyama method 
[9]. The algorithm may be described as follows: 

Algorithm l. Let r and cr be the constants of the problem, a I the time-step, Tol ihe 
tolerance, and nmn / he number of simulations* Then, we have: 

for each m in I : nmn 

Sian with ihe initial value Sq 
for each n in 0, 1,2... 
if $„ < Tol 

extiiiction_time(m) *— nAl 
break 
end if 

Generate a random number w with normal distribution fsf( 0, 1) 

^ S„ + pS„At + crS n \fAtw 

end for 
end for 

Compute the mean and standard deviation of ext incti on_time 
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In Figure 1, we have plotted in black the function t( ~s) = ln( ~s + 1)/ 1 [3 1 , for p = 0.01,a = 1. The red circles located on 
the black curve correspond to the results from Algorithm 1, executed with time step At = 10-3, tolerance Tol = 1, 
nrun = 105 and 

S (0) = 10, 50, 100, 200, 300, 400. Observe that we have given Tol and the same values, i.e., Tol = = 1, to make 

both approaches comparable; indeed, their matching is obvious from the figure. 



FIGURE 1 COMPARISON BETWEEN r (T) = In (s + 1)/\J3\ (BLACK CURVE) VERSUS THE RESULTS FROM ALGORITHM 1 (RED 
CIRCLES). THE MATCHING OF BOTH APPROACHES IS REMARKABLE. 


Conclusions 

The results obtained with both techniques are largely coincident, it is very important to underline that each 
execution of Algorithm 1 may take run time (and it has to be repeated again for each initial value), while function 
t( ~s) is explicitly known for any point. 

Moreover, when studying a real market price the volatility is not constant as assumedin the Black-Scholes model, 
but varies. In this way, many authors have proposed that the volatilities should be modeled by a stochastic process 
[6]. As already commented in [4], the mean extinction time can also be computed by solving the stationary 
backward Kolmogorov equations [2]. The equations are linear second-orderpartial differential equations with 
variable coefficients and, therefore, it is impossible in general to find exact solutions, so we can only compute 
numerical approximations. 

On the other hand, recent results in [10] and [5] about convergence rates for the approximations 

these stochastic volatility models, in my opinion, will allow to estimate and compare mean extinction-time in these 
more complicated stochastic models. 
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